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NOMENCLATURE 


D 

E. 

H 

J 

L 

S 

T 

V 
W 
W 
X 

Y 
Z 


Orifice diameter 

Energy exchange efficiency 

Duct height 

Momentum flux ratio, (pV^)j/(pV^) w 

Vertical distance from jet centerline 

Orifice spacing 

Temperature 

Velocity 

Width 

Weight flow rate 

X distance, axis parallel to Primary 
Y distance, vertical axis 
Z distance, lateral axis 


Subscripts 

c 

cent 

EB 

i 

j 

mid 

min 

v 


1/2 


Temperature centerline 
Jet centerplane (Z - O.O) 
Energy balance value 
Properties at a point 
Jet condition 
Jet midplane (Z = S/2. ) 
Minimum value 
Velocity centerline 
Primary stream condition 
Half value or half width 


Superscripts 

+ Plus side of jet (far side from injection plane 

Minus side of jet (side near injection plane) 



Nomenclature (cont.) 


Greek 

0\ Temperature difference ratio, (T^ - T^)/(T W - Tj) 

Energy balance temperature difference ratio (T^ - T£g)/(T M " ^j) 
p Density 

® ^c.cent’ 


viii 



I SUMMARY 

The purpose of this study was to correlate the experimental diluent 
ai r/primary combustor gas mixing efficiency and downstream temperature dis- 
tributions obtained during the Multiple Jet Study (NAS 3-15703) to gas turbine 
combustor operating and design variables. The experimental data were generated 
by probe measurements from tests on single rows of multiple dilution orifices 
(diameters of .64 to 2.54 cm) injected into a low Mach number (M = .03) heated 
primary stream (450°K to 750°K) in a 10.2 by 30.5 cm duct. The correlations 
were developed using power form or exponential equations which related the 
various dependent temperature field variables to the independent operating and 
design variables. 

The dependent mixing and jet penetration parameters correlated at 
each downstream data location were: the jet/primary stream mixing efficiency; 
the jet temperature and velocity trajectories downstream of the injection ori- 
fice in the jet centerpl ane-of- symmetry; the maximum centerplane temperature 
difference (which is on the temperature centerline); the jet half-width values 
on each side of the jet centerline in the jet centerplane-of-symmetry; and 
the minimum temperature difference values on each side of the centerline. 

When coupled with the Gaussian form assumed for the profiles, these parameters 
completely define the centerplane temperature distribution at any downstream 
location. 

The development of the off-centerplane temperature distribution made 
use of the observed Gaussian nature of the vertical temperature distribution 
at all stations where the flow field was influenced by the diluent jets. The 
off-centerplane correlations included the ratio of temperature maximum values 
at the lateral planes to the maximum values in the centerplane and the ratios 
of jet thermal penetration in the lateral planes to the thermal penetration 
in the centerplane. The off-centerplane half-widths were assumed to be equal 
to the corresponding centerplane values. Also the ratio of the minimum to 
maximum temperature difference at any off-centerplane location was assumed 
to equal the corresponding centerplane ratio. 
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I Summary (cont.) 

The parameters and relationships described above provided the necessary 
input for complete characterization of the temperature field downstream of the 
diluent injection plane. The range of the operating and design variables used 
to develop the various correlations were selected to make the correlations 
relevant for use in the design of a wide spectrum of combustors for gas tur- 
bine engines. Five independent variables (one operating variable and four 
geometric variables) were used to achieve correlation of the dependent para- 
meters with the test data. The independent variables used in the correlations 
and their ranges were: jet/primary momentum flux ratio, (pV )j et /(pV ^ primary 
(6-60); orifice spacing/jet diameter ratio, S/D j (2. 5-7. 5); duct height to 
jet diameter ratio* 'H/Dj (5-20); downstream distance to jet diameter ratio, 
X/Dj (1.25 - 30); and lateral distance to jet spacing, Z/S (0-.5). In addi- 
tion, diluent to primary flow ratios of .04 to .60 were implicit in the data 
but were not required to correlate the data. The correlations were based on 
data obtained from a matrix of five axial stations, six lateral stations and 
20 vertical stations in the flow field during approximately 50 tests on eleven 
orifice row designs. 

II INTRODUCTION 


The “Program to Correlate Diluent Air/Primary Combustion Gas Mixing 
Parameters with Gas Turbine Operating and Design Variables'*, was conducted 
under NASA Lewis Research Center contract NAS 3-18026. The correlations 
developed were based on data generated during the Multiple Jet Study (Contract 
NAS 3-15703, Ref, 1). A mixing efficiency parameter, termed the energy exchange 
efficiency (E^)* was defined during the performance of the Multiple Jet Study 
and was shown to quantify the diluent/primary stream mixing efficiencies over 
a range of test and operating conditions (Ref. {!)). This study extended the 
previous study to mathematically define the relationship between and the 
combustor variables. Also, the study included an investigation of the corre- 
lation between the combustor variables and the temperature profiles downstream 
of the diluent injection plane. A goal of the program was to develop a general 
model which would allow predictions of flow field temperature distributions 
as a function of combustor operating and design variables. 
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II 


Introduction (cont.) 


Correlations of the penetration and mixing of jets in a crossflow 
has application to many problems of current interest, such as: 

(1) Cooling of primary combustion qases with diluent air in 
gas turbine combustors. 

(2) Cooling of hot gas streams in numerous industrial and military 
devices. 

(3) Film Cooling of combustion chamber walls, turbine blades, 
and reentry vehicle nose cones. 

(4) The aerodynamics of STOL and VTOL aircraft. 

(5) The concentration and paths of pollutants downstream of 
industrial chimneys or downstream from discharge lines 
leading into rivers or streams.- 

The results of this study apply most directly to Items (1) and (2) 
above. The development of valid correlations for the mixing process between 
cool multiple jets and a hot primary gas stream has two principal interrelated 
benefits: (1) through proper design of secondary air admission ports, the 
combustor weight is reduced and packaging is improved since lengths required 
to achieve uniform temperature and mass flux profiles can be minimized, and 
(2) the decreased combustor length required for complete mixing will result 
in minimum residence time for production of nitrogen oxides. 

Although the interaction of subsonic circular and noncircular jets 
injected normally into a subsonic mainstream flow has been the subject of 
numerous analytical and experimental studies, (Ref. 2-7), most published 
works to date have dealt with single jets rather than multiple jets in a 
bounded cross flow as required to simulate the gas turbine combustor secondary 
air admission problem. Two recent experimental studies, the above mentioned 
work done by Aerojet Liquid Rocket Company (Ref. 1) and work done by Case 
Western Reserve University (Ref. 8), have produced data for the study of the 
interaction of a row of multiple jets in a confined crossflow. Correlation 
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II 


Introduction (cont.) 


of a portion of the Reference 1 data has been done by Cox at Pratt & Whitney 
Aircraft (Ref. 9 and 10). The present study is based on a larger body of 
data than the Reference 9 study and the correlations were derived over a 
wider range of variables. 

The correlations presented here were developed by relating the 
various dependent temperature field variables to the independent operating 
and design variables using power form or exponential equations. The basic 
forms of the correlating equations were developed from theoretical considera- 
tions and from observations of the empirical behavior, with the specific 
coefficients and exponents derived from a covariance analysis of the test 
data. This technique has led to correlations which are simple to apply 
and lead to an insight into the physical processes occurring during pene- 
tration and mixing of multiple jets in a confined crossflow. 

Ill TECHNICAL DISCUSSION 

A. Data Sample and Method of Analysis 

The multiple jet correlations are based on data obtained during 
the performance of Contract NAS 3-15703 (Ref. 1). The centerplane correla- 
tion equations are based on multiple covariance analyses using over 200 test 
data points from eleven orifice row configurations at an average of four test 
operating conditions. A summary of the test configurations and operating 
conditions is contained in Table I. For the off centerplane evaluation, data 
from over 800 test data points was used. Although the correlations were based 
on the Reference 1 data, some comparisons are made with the experimental results 
of Reference 8. In addition, the results of the present study are compared 
with the results of Reference 9, which was based on selected tests from the 
Reference 1 data. 

The covariance analyses were conducted using ALRC One-Way Multiple 
Covariance Analysis Program (FD 0088). The program uses standard multiple 
regression and covariance techniques and computational methods. The analysis 
may be performed for up to 20 variables and 500 groups. A trans-generation 
feature allows for additional variables to be generated or transformed from 
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Ill Technical Discussion (cont.) 


the input variables as desired. For the particular requirements of the 
"Multiple Jet Correlation Study" the program was modified to accept input 
from the mass storage data files created for each dependent parameter. 

Five independent variables (one operating variable and four 
geometric variables) were used to achieve correlation of the dependent para- 
meters with the test data. The independent variables were: jet/primary 
momentum flux ratio, (p V 2 )j et /(pV 2 ) iniary ; orifice spacing/jet diameter ratio, 
S/D j; duct height to jet diameter ratio, H/Dj ; downstream distance to jet 
diameter ratio, X/D,; and lateral distance to jet spacing, Z/S. The diluent 
jet to primary stream density ratio was an additional parameter which was 
varied during the test series. However, over the range of density ratios 
tested (1.6, 2.2 and 2.7), no significant influences of the parameter were 
observed. Not used as a parameter to correlate the data, but implicit in 
the data, were diluent to primary flow ratios of .04 to .60. The correla- 
tions were based on data obtained from a matrix of five axial stations, six 
lateral stations and 20 vertical stations in the temperature field. The 
ranges of the operating and design variables used to develop the various 
correlations are given in Table I. 

B. Correlating Parameters and Assumptions 

1. Mixing Efficiency 

A mixing efficiency parameter, termed the energy exchange 
efficiency (Ej) was defined during the performance of the Multiple Jet Study 
(Ref. 1) and was shown to quantify the diluent/primary stream mixing efficiency 
over a range of test and operating conditions (Ref. 1, 11). During the present 
study the Ey values were correlated as a function of the downstream distance, 
the combustor momentum flux ratio and the diluent jet size and spacing. The 
advantage of developing a correlation for a general mixing efficiency para- 
meter, such as Ey, in addition to the other temperature field parameters, is 
that evaluation of this single parameter will allow the designer to quickly 


*The orifice spacing to duct height ratio, S/H, also proved to be a valuable 
independent parameter, and was used in place of S/D^ in two of the correlations. 
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Ill Technical Discussion (cont.) 


estimate the overall efficiency of the diluent/primary stream mixing process 
without the need to evaluate the many separate equations necessary for complete 
temperature field analysis. 

2 . Temperature Field Parameters 

An illustration of the coordinate system used during the 
study together with a representation of the temperature field parameters is 
contained on Figure 1. In order to define the dimensionless temperature 
field downstream of the diluent injection orifices correlations must be 
developed for certain principal parameters and some key observations regarding 
the nature of the temperature field must be utilized. To develop the temp- 
erature field the diluent jet temperature trajectory in the orifice center- 
pi ane-of- symmetry downstream of the injection orifice must be defined and the 
temperature values along this path must be known. In addition, vertical 
temperature distribution shape parameters in the orifice centerplane must 
be defined (See Figure 1) and the shape of the temperature distribution off 
the centerplane must be known. 

(a) Jet Trajectory Parameters 

The diluent jet trajectory is defined in terms 
of the local penetration depth as a function of downstream distance, with 
both the penetration and downstream distance nondimensional ized by jet dia- 
meter. Both a velocity penetration and a thermal penetration were evaluated 
during this study. The velocity penetration, Y^/Dj, is defined as the loca- 
tion of the maximum total pressure. The thermal penetration, Y^/D^, 1S defined 
as the location of the maximum temperature difference. The locus of penetration 
with downstream distance defines the trajectories. The thermal penetration has 
a direct impact on subsequent correlations for the complete temperature field. 

(b) Non-Dimensional Temperature Parameters 

The temperature parameter used for this study is 
the nondimensional temperature difference in the flow field downstream of jet 
injection, theta (O'-), defined as: 
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Ill Technical Discussion (cont.) 



( 1 ) 


where: 

O'. = Theta, nondimensional temperature difference at a 

1 point in the flow field 

T = primary flow stagnation temperature 


T. 


jet stagnation temperature 

stagnation temperature at a point in the flow field 


Theta is a measure of the temperature suppression in the flow field. The 
value of theta can vary from one, when measured temperature equals the jet 
temperature, to zero, when the measured temperature equals the primary stream 
temperature. 


If complete mixing of the jet and mainstream 
flow occurs, the value of theta will be constant and T^ will be everywhere 
equal to the ideal equilibrium temperature between jet and mainstream. Thus, 




where: 


( 2 ) 



ideal equilibrium theta 

stagnation temperature resulting from complete 
thermal energy exchange 


The ideal theta is a useful parameter; a comparison between the measured 
local theta and the ideal theta provides a means of gauging the local mixing. 


The maximum dimensionless temperature difference 

on the centerplane, ft, cent » defines the thermal trajectory. For the case of 

a single jet in a semi-infinite crossflow, cent and ® c ,cent is 

expressable as O' . Ref. 7. For multiple jets in a confined flow, 

c , cent 
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1 ^c.cent^ %’ and the power form is not appropriate. If the centerline 
temperature decay is expressed as, 




T c,cent 

- T. 

J 


(3) 


® is a measure of the flow field temperature reduction occurring along the 
centerline compared to the maximum possible reduction. Since 1 s(h)sO s 
® can be modeled with the power form. Then ^ cent can be obtained from 


& 


c,cent 


= ® 0 - fr EB ) + cr ( 


EB 


(4) 


(c) Centerplane Temperature Profiles 


The correlation of the thermal trajectory and the 
centerline maximum temperature difference are the first steps in a system of 
equations to define the flow field temperature distributions. The next step 
is the determination of the temperature profile shape factors which will allow 
the temperature distribution in the orifice centerplane about the jet centerline 
to be predicted. From the work of Holdeman (Ref. 7) and Cox (Ref. 9) and the 
data of Reference (1), the assumption of a Gaussian vertical temperature dis- 
tribution appears to offer a simple yet accurate means of modeling the data. 


Here another important difference between the 
single jet flow and the multiple, confined jet flow must be recognized. That 
is, 0 does not have to decay to zero with increasing radial distance from the 
centerline. Thus the minimum dimensionless temperature difference above and 

4 - 

below the centerline, fr— . may be greater than zero, and must be 

mi n , cen l 

correlated. Also, the traditional definition of the half width (the width 
where fr = O' ,/2) must be modified such that W — , /9 /D. is the distance 

C j C6n t l/c J 

from the centerline to where O' = (0 cent + Cen t^ 2 * is necessar y 

since O' may be everywhere greater than cent /2, and the traditional half- 
width would be undefined. 


Using these parameters, the vertical temperature 
distribution in the centerplane was defined by: 
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Ill Technical Discussion (cont.) 


r - 



where: 


AO. 

(V 

i ,cent 

-0 

+ 

"min, cent 

II 

<l 

fr' 

c,cent 

-Cr 

+ 

min, cent 


L./D, = local distance from centerline nondimens ional ized by 

1 jet diameter 


W 



plus or minus side half width nondimensionalized by 
jet diameter 


A schematic drawing of the test duct is shown on Figure 1 with a typical 
vertical centerplane temperature profile and temperature field parameters 
illustrated. 


(d) Lateral Plane Temperature Profiles 

Correlations for the vertical temperature dis- 
tributions off the centerplane were needed in order to model the complete 
temperature field. These off-centerplane correlations included the ratio 
of the maximum temperature difference values at the lateral planes to the 
centerline values in the centerplane (Ch. z /0" c ^.) and the ratios of the jet 

thermal penetration in the lateral planes to the thermal penetration in the 
centerplane (V /Y .). In addition to these correlations, the development 
of the off-centerplane temperature distributions made use of the observed 
Gaussian nature of the vertical temperature distribution at all stations where 
the flow field was influenced by the diluent jets. Also, the observation that 
the ratios of the minimum to maximum temperature difference at any off-center- 
plane location were essentially equal to the corresponding centerplane ratios 
was a key modeling relationship used in defining the complete temperature field. 
Another major simplifying assumption, justified by the experimental data, was 
that the off-centerplane half-widths were equal to the corresponding center- 
plane half-widths. 
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Ill Technical Discussion (cont.) 


The parameters and relationships described in 
the preceding paragraphs provided the necessary input for complete charac- 
terization of the temperature field downstream of the diluent injection plane. 
The correlations were developed by relating the various dependent temperature 
field variables to the independent operating and design variables using power 
form or exponential equations. The basic forms of the correlating equations 
were developed from theoretical considerations and from observations of the 
empirical behavior, with the specific coefficients and exponents derived 
from a covariance analysis of the test data. A summary of the correlation 
equations is shown in Table II. 

C. Mixing and Centerplane Correlation Equations 
1. Energy Exchange Efficiency 

An energy exchange efficiency parameter was defined 

in Reference 1 by: 


E 


T 



r 



100 


w 


T 


where: 


Wj.| - local jet mass flow rate 

W o3 .j = local primary mass flow rate 


W 


T 


total mass flow rate 
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Ill 


Technical Discussion (cont.) 


The energy exchange parameter expresses the mixing effectiveness (in percent) 
as the energy exchanged between the cool jets and the hot primary stream, at 
any axial station, compared to the energy exchanged if both streams came to 
thermal equilibrium.. The E^. values have been shown to quantify the diluent/ 
primary stream mixing efficiencies over a range of test and operating condi- 
tions (Ref. 1, 11). During this study the energy exchange parameter has been 
correlated to the combustor operating and design variables by the following 
relationship: 

-a 

E-j. = TOO [1 - e ] ( 6 A) 

where: .41 .44 -1.0 .44 

a = .682 J (S/Dj) (H/Dj) (X/Dj) ■ ( 6 B) 

A plot of the E t correlation equation, which has a one 
sigma standard error of prediction of 5.6 is shown on Figure 2. Inspection 
of Equations ( 6 A) and ( 6 B) shows Ej to be bounded by values of 0 and 100 and 
shows the E-j. prediction to increase with increasing momentum flux ratio, J, 
orifice spacing, S/Dj, and downstream distance X/Dj* and orifice size 1/(H/Dj). 
The correlation was developed over the ranges of independent variables given 
in Table I, but excluded those specific cases (approximately 10 % of the data) 
where jet over penetration occurred, i.e., cases combining high momentum flux 
ratio with large hole size and hole spacing. 


2. Jet Velocity Penetration 


V D j> was: 


Y /D 

V 


J 


The correlation obtained for the jet velocity penetration, 


.12 .23 .57 .18 

.549 J (S/Dj) ( H/D ) j (X/Dj) 


(7) 


From the form of Equation 7 one may see that increasing momentum flux ratio, 
duct height/orifice diameter and/or spacing increases the trajectory path 
depth. The agreement between the data and the correlation is shown on Figure 3. 
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Ill Technical Discussion (cont.) 

Approximately 86% of the data are within a + 20% band about the prediction line. 
This data band is a consequence of the very uniform vertical velocity distri- 
bution shown by a large portion of the data. The uniform velocity distribution 
caused some random scatter in the location of the maximum velocity values, how- 
ever, the covariance analysis indicated good correlation with all of the above 
independent variables. 

Velocity penetration data was also available from the 
work of Kamotani and Greber (Ref. 8). These data indicate less jet velocity 
centerline penetration than is predicted by Equation (7), except at the highest 
momentum flux ratios when the data from Reference 8 shows greater penetration 
than does the prediction. The data from Reference 8, for tests with H/D = 8 
and S/D = 2, is shown on Figure 4, along with the corresponding trajectory 
predictions using Equation 7. Differences in primary stream boundary layer 
effects and jet velocity profiles may partially account for these penetration 
differences shown on the figure. The jet velocity profiles from Reference 8 
corresponded to fully developed pipe flow while the Reference 1 work used sharp- 
edged orifices and the jet velocity profiles were not fully developed. Jet 
velocity profile differences between pipe flow and nozzle (or orifice) flow 
were observed to cause approximately a 10% reduction in jet penetration for the 
pipe compared to the nozzle (Ref. 8). If the corrections for velocity profile 
and boundary layer development are made to the predictions on Figure 4 agree- 
ment between measured and predicted values is improved at the lower momentum 
flux levels, but is worse at J = 72. The variation of the trajectory with 
downstream distance appears to be correctly given by Equation (7). 

For most of the data surveyed the agreement between the 
Reference 8 data and the predictions of Equation (7) appeared best at a momentum 
flux ratio of 32. For much of the Reference 8 data low momentum flux ratios 
(J = 8) resulted in substantially less penetration than did the data of Reference 
1, upon which Equation 7 is based. At high J values the Reference 8 data shows 
more penetration than does that of Reference 1. Apparently the influence of 
momentum flux ratio on jet penetration from the two sets of data are significantly 
different. A log-log plot of the penetration depth as a function of momentum 
flux ratio is shown on Figure 5 for both the Reference 8 data and the Reference 1 
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Ill Technical Discussion (cont.) 


data with two orifice row configurations* S/D = 2 at H/D = 8 and H/D = 12. 

The data is shown at a location 10 diameters downstream of the injection plane, 
The data from Reference 1 have a constant exponent on 0 while the Reference 8 
data indicate an increasing exponent on J with increasing J. 


3. Jet Thermal Trajectory 


Y c /D JS was: 


where: 

b = 


The correlation obtained for the jet thermal penetration, 


.25 .14 .38 .17 -b 

Y c /Dj = .539 J (S/Dj) (H/Dj) (X/Dj) e (8) 


(X/H) 2 (H/S - \fj/3.5)/11.0 


(9) 


As with the velocity trajectory, increasing momentum flux ratio, duct height/ 
orifice diameter and/or orifice spacing all tend to increase the depth of the 
trajectory path. However, for the thermal trajectory an exponential modifier 
is used to model path recurving which occurs with under penetration at far 
downstream distances. A correlation for Y /D , was derived by Cox in Reference 

C u 

9. The Reference 9 correlation is based on a baseline data case with corrections 
to the baseline case obtained from polynomial (up to 4th order) curve fits on 
X/D,. Comparison of the correlation equation (8) with the Reference 9 correla- 
tions showed Equation 8 matched the data slightly better than do the Reference 
9 correlations. The correlations of Reference 9, due to the polynomial curve 
fits, are not applicable for X/D ^ 21. 

■J 

The agreement between the data and correlation of Equation 
8 is shown on Figure 6. As with the velocity trajectory, the thermal trajectory 
definition was difficult due to the uniform vertical temperature profiles of a large 
portion of the data. Approximately 85% of the data falls in a + 20% band about 
the prediction. At the far downstream locations the data scatter is more evident 
than at locations near the orifice injection plane. The covariance analyses 
indicated significant exponents for all the specified independent variables. The 
validity of the trajectory equation is evidenced by the good agreement between 
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Technical Discussion (cont.) 


measured and predicted temperature profiles which will be shown in Section HIE'. 
4. Jet Centerline Temperature Difference Values 


The correlation obtained for the jet temperature centerline 


values was: 


(V 


c,cent 


-.4- 


1.536 J 


X/D 


1.15 


(i - o^g) + cr EB 


(io) 


where: 


EB 

f 


the ideal theta defined in Equation 2 




S/H / (1 + S/H) 


From Equation 10 the temperature centerline values, 0^ decrease with 
increasing downstream distance and momentum flux ratio and is strongly 
influenced by 0^. Also the influences of X/D j and J on cgnt are coupled 
to the spacing, S/H. The agreement between the measured data and the corre- 
lation Equation 10 is shown on Figure 7. The data on Figure 7 are shown plotted 
as the prediction value as a function of the measured value, since a single 
correlation curve as a function of X/Dj can not be drawn due to the variable 
power on X/Dj in Equation 10. Approximately 85% of the data falls in a + 

10% band about the correlation line. 


Centerline temperature difference ratios were measured 
for heated jets injected into a cool primary stream in the work done by Kamotani 
and Greber in Reference 8. The rates of change of the dimensionless temperature 
ratio, ©^ c e n t ^ as a f unc ^ on downstream distance for the Reference 8 data 
were approximately the same as that shown by the cool jets in heated primary 
stream data used on this program. A correlation for the jet centerline dimension- 
less temperature ratio based on a portion of the Reference 1 data was presented in 
Reference 9 as an exponential decay. The form of the Reference 9 equation differed 
from the more conventional power form, and the prediction appeared to diverge 
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Ill Technical Discussion (cont.) 


from the measured data at large X/D J5 although the limits on 0" c cent were 
well defined. 


5. Plus and Minus Side Minimum Temperature Difference Values 


As mentioned previously, recent studies (Ref. 1, 7, and 9) 

have shown the vertical temperature distribution in the orifice centerplane 

to be approximately Gaussian in nature. Therefore the distribution can be 

modeled if the location (Y cen ^./Dj) and magnitude (0^ cen -j.) .of the peak theta 

values are known and if the distance from the centerline to some characteristi c 

theta values (such as a half value) on the near (-) and far ( + ) injection sides 

of the jet centerline can be defined. For the case of single jet injection 

the characteristic distance dimension is from the centerline to the theta half 

values, O' J2. For multiple jet injection temperature difference as low 
c,cent + 

as # c cent /2 may not exist on the centerplane. Thus the half-widths, w-i /2 /d., 
are defined as the distance from the centerline to the location where: 




1/2, cent 


^®c,cpnt + ^min,cent^ 


OD 


To specify the profile using this definition of the half-width, the ratio, 

(0^- . J/CCT J must be known for all conditions. The form chosen for 
' mm, cent" c,cent 

these correlations was: 

+ 


+ 

1/2, centre, cent 



(12a) 


+ 

for O' 


1/2, cent* 


and the corresponding form for the minimum value: 

+ 


+ -c 

ty . ./O' , = 1 -e 

min, cent c,cent 


(12b) 


For the plus side ratio: 

+ - 1.62 1.5 -3.67 1.1 

C = 0.038 J (S/D ) (H/Dj) (X/Dj) (13) 
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Ill Technical Discussion (cont.) 

This correlation results in increasing plus side minimum theta ratio values 
with increasing downstream distance, X/D J5 increasing momentum flux ratio, J 9 
and increasing jet spacing, (S/Dj), and jet diameter, (H/Dj)" . These results 
are reasonable because increasing all the above mentioned parameters would 
increase jet penetration and thus result in a trend toward higher plus side 
theta minimum values. 

For the minimum theta values on the minus side of the 
jet the correlating function, c 5 in Equation 12 was: 

-.3 -1.4 .9 

c = 1.57 J (S/Dj) (X/Dj) (14) 

This correlation predicts increasing minus side theta ratios with increasing 
downstream distance, but with decreasing momentum flux ratio and orifice 
spacing. The orifice size did not significantly influence the minus side 
minimum theta ratio. The inverse relationship between the changes in the 
minimum theta ratio and changes in momentum flux ratio and spacing is probably 
due to the fact that jet penetration increases with J and S/D js which would 
allow the jet minus side theta values to decay to lower minimum values. The 
agreement between the data and the predictions for the plus and minus fly 2 , cent 
values are shown on Figures 8 and 9 respectively. 

6. Plus and Minus Side Half Widths 

+_ 

As discussed in the preceding paragraph, the O' y2,cent va ^ ueS5 
(Eq. 11), were the dimensionless temperature parameters used to define a charac- 
teristic dimension, the half width, used in the Gaussian dimensionless tempera- 
ture distribution equation (Eq. 5). The correlation for the plus side half width 
nondimensionali zed by the jet diameter, Dj, was: 

+ .18 -.25 .5 .5 

“ 1.2,cent /D 0 = - 16Z J (S/D J } {H/D J ) (X/D 0 ) (15) 

The correlation equation for the minus side half width was 
derived by difference from correlations of jet total half width and plus side 
half width. The resulting correlation was: 
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Ill Technical Discussion (cont.) 

W ’l/2,cent^J = -2T 15 (S/Dj )' 27 H/d/ (X/Dj)' (16) 

Difficulties encountered in a direct correlation of the minus side half 
width were probably a consequence of the very uniform minus side dimension- 
less temperature profiles for a large portion of the data. This made defini- 
tion of the precise location of the minus side minimum theta values difficult. 

The half width correlations can not, by themselves, be 
related to changes in the dimensionless temperature profiles since the half 
width values must be coupled with the corresponding minimum and centerline 
theta values in order to properly interpret the influences on the dimensionless 

temperature profiles. For example, if ®^ m -j n5Cen t and ®c,cent are near ^ e q ua l 
a uniform temperature profile will result, even for very small half width 

values. 


D. Off-Centerplane Correlations (Z Planes) 

Two off-centerplane correlation equations were developed: (1) the 
ratio of the maximum temperature difference at the lateral (Z) planes to the 
centerline values in the centerplane and; (2) the ratios of the jet thermal 
penetration in the lateral (I) planes to the thermal penetration in the center- 
plane. The observed Gaussian nature of the vertical temperature distribution, 
at all stations where the flow field was influenced by the diluent jets, was 
used to define temperature field profiles at the off centerplane locations. 

The data showed the ratio of theta minimum to theta centerline values at any 
location off the centerplane were essentially equal to the corresponding center- 
plane ratios. Thus the previously developed centerplane minimum theta correla- 
tions could be applied at the off centerplane locations. Also, the off center- 
plane half-widths were assumed to be equal to the corresponding centerplane 
hal f-widths. 


1. Ratios of Maximum Theta Values in Lateral Planes 

to Theta Centerline Values in Centerplane, & c z /& CyCen t 

The basic form of the correlating equation for the lateral 
plane to the centerplane theta ratio z /°c,cent Was: 
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where: 

Z 

S/2 


r /O' 
c,z' c,cent 


] 


& 

c,mid _ 
O' 

c ,cent 



local distance from centerplane to plane Z 
distance from centerplane to midplane 


( 17 ) 


This form makes use of the mid to centerplane theta ratios and the lateral 

position ratio, Z/(S/2). Using Equation 17 the predicted theta ratios will 

be between 0 and 1 and the rate of change of & with 1 will go to zero at 

c , z 

the centerplane. The power on Z will cause the variation of the theta ratio 
with lateral distance to be parabolic. A better basic form might be one which 
will allow the variation of the theta ratio with lateral distance to contain 
an inflection point and have zero slopes at both the centerplane and midplane. 
However, at the present time this more sophisticated modeling doesn't appear 
justified or necessary. The correlation equation for ^^-j^^c^ent is: 


-d 


) /& = 1 - 
c,mid / c,cent 


(18) 


where: 


d 


.53 -1.53 .83 

.452 J (S/Dj) (X/Dj) 


(19) 


Thus the midplane to centerplane ratio increases with increasing momentum 
flux ratio and downstream distance (more jet spreading) and decreases with 
increasing orifice spacing. The dimensionless jet diameter, [H/Dj] ^ did 
not appear to significantly influence the theta ratios. 


2. Ratio of Penetration Depth in Lateral Planes to 
Penetration Depth on Centerplane, Y c Z /Y c CGnt 

The basic fornj of the Y c 2 /Y c cent correlation was 
identical to that used for O’ /© t ; 

L 5 4 C 5 L c 1 1 L 
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with 


Y /Y 
c,z c,cent 


r 



V 

c ,mi d 
y 

c,cent 



( 20 ) 


Y /Y 
c ,mi c,cent 


" 9 

e 


( 21 ) 


where: 

9 


.67 -1.0 .54 

.227 J (S/Dj) (X/Dj) 


( 22 ) 


The trends predicted by this correlation are similar to those predicted by 
the theta ratio correlation equations. 

E. The Complete Temperature Field 


The parameters and relationships described above provide the 
necessary input for complete characterization of the temperature field down- 
stream of the diluent injection plane. A computer code, FIELD, was developed 
which incorporated the various equations and relationships into a temperature 
field model. These correlation equations were summarized in Table II. A listing 
of' this code is contained in the Appendix along with a sample input. The predicted 
temperature profiles for Figures 10 through 29 were obtained using the FIELD 
program. 


1. Dimensionless Temperature Profiles in the Centerplane 

Predicted and measured dimensionless temperature profiles 
in the orifice centerplanes are shown on Figures 10 through 20. The test 
configuration matrix of orifice sizes and spacings used to develop the center- 
plane correlation equations are shown on Table III along with the momentum flux 
ratios surveyed. The specific configurations and momentum flux ratios selected 
for centerplane profile illustration are shown on Table IV. The centerplane 
profiles are shown for downstream distance to duct height ratios of .25, .50, 
1.0 and 2.0, with the exception of the H/Dj = 20 case where X/H values of .125, 
.25, .5 and 1.0 were used. 
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Figure 10 contains data from a H/Dj = 10.2 and 5/D j = 
3.8 orifice row configuration at a momentum flux ratio of 26.7. This con- 
figuration approximates an ‘'average" configuration based on H/Dj and S/Dj 
values. Good agreement between the empirical data and. the pediction may be 
seen at all four downstream planes. 


Data obtained with H/D = 10 and momentum flux ratio 
of 26 is shown on Figures 11, 12 and 13 for orifice spacings, S/Dj, of 2.5, 
5.1 and 7.7 respectively. These data show the predicted and measured in- 
creases in jet penetration as spacing is increased. Agreement is again good 
between the experimental data and the prediction except for the S/Dj =7.7 
case at the first two planes when the penetration depth is under predicted 
by approximately 10%. The data in Figures 11 and 12 were used in demon- 
strating the correlation method of Reference 9, and the predictions from 
this reference are shown for comparison. 


The effect of momentum flux ratio on the predicted and 
measured dimensionless temperature profiles are shown by the data of Figures 
14 and 15 for nominal H/Dj « 10 and S/Dj =5 at nominal momentum flux ratios 
of 6 and 62 respectively. The data of Figures 14 and 15, along with the 
J = 26 data of Figure 12 show the increase of jet penetration with momentum 
flux ratio. 


The data from tests of the smallest orifices, H/Dj = 20, 
at the closest spacing, S/Dj = 2,5 is shown on Figure 16, for a nominal 
momentum flux ratio of 25. Both the measured and predicted data show the 
small penetration distances achieved at all stations. Agreement between the 
prediction and the measured data is very good at the three upstream stations 
but only fair at X/H = 1.0. Figures 17 and 18 contain data from tests using 
a nominal H/D^ of 15 at S/Dj values of 2.5 and 5.1, respectively, and nominal 
momentum flux ratios of 60. Again agreement between the empirical data and 
the predictions appears good at most stations and the increase in S/Dj is 
shown to increase jet penetration. 
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Ill Technical Discussion (cont.) 

Comparison of the measured and predicted profiles for 
the larqest orifice diameter tested (H/D. - 5) is shown on Figures 19 and 20 

' u 

for J = 13.3 and 5/Dj =2.5 and for 0 = 27.2 and S/Dj = 5, respectively. 

The prediction for the 0 - 13.3 test appears to match measured data well. 

For the 0 - 27.2 case the prediction for X/Dj =1.3 and X/Dj = 2.5 under- 
estimates the jet penetration. 

The test conditions used to illustrate the model applica- 
bility on Figures 11, 12, 15 and 19 were also used in the study of Reference 
9. The Reference 9 predictions are shown on the figures for comparison. 

Based on these data the centerplane predictions using the correlations from 
this study appear to model the empirical data as well or better than do the 
predictions of Reference 9. In addition, the simplicity of the correlations 
developed during this study allows easy computation, provides some insight 
as to the physical processes occurring during penetration and mixing, and will 
allow confident extrapolations. 

2. Dimensionless Temperature Profiles in the Lateral (Z) 

PI anes 

Predicted and measured dimensionless temperature profiles 
in the lateral planes are shown on Figures 21 through 29. The test configura- 
tion matrix of orifice sizes and spacings used to illustrate the lateral 
plane profiles are shown in Table V. The lateral planes shown on the figures 
are for Z/S = 0.0, (centerplane) , Z/S = .2, Z/S ~ .3, and Z/S = .5 (midplane). 
With the exception of Figure 28 which shows data at X/H = .25 all the profiles 
are shown at a downstream distance to duct height ratio of 1.0. The data shown 
on Figures 21, 22 and 23 are for nominal H/D j = 10, S/Dj = 5 and nominal 
momentum flux ratios of 6, 27 and 62 respectively. Both the measured and pre- 
dicted data show the increase in jet penetration and the increasing spreading 
of the jet (less profile change with Z lateral plane) as momentum flux ratio 
is increased. 
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Data from tests with a nominal H/D^ = 10 and nominal 
momentum flux ratio of 26.0 are shown on Figures 24 and 25 for S/D, values 

u 

of 2.5 and 7.7 respectively. Comparison of these data and the data of 
Figure 22 shows the increase in centerplane jet penetration, and the flattening 
of the temperature profiles in the Z lateral planes, with increasinq'orifice 
spacing. Good agreement between predicted and measured temperature profiles 
is evident on Figures 21 through 25. 

Lateral plane temperature profiles for the smallest 
jet diameter, H/Dj = 20 and smallest spacing S/D^ =2.5 are shown on Figure 
26. These data show the flat minus side temperature distribution in both 
the Y and Z directions. For the plane shown, X/H = 1.0, the predicted pro- 
files underestimate the jet penetration slightly, agreement is better at the 
upstream stations as may be seen from the centerplane data of Figure 16. 

Predictions for operating and design conditions used 

in the study of Reference 9 are shown on Figures 27, 28 and 29, with the 

predictions from this reference (or Reference 10) shown for comparison. 

The data of Figure 27 are for X/H = 1 , J = 57.3, H/Dj = 15 and S/Dj = 2.5. 

Figure 28 shows data at X/H = .25, J = 24.7, H/Dj = 15 and S/Dj = 5 and the 

Figure 29 data are for the largest orifice tested, H/D, = 5, at X/H = 1.0, 

j 

5/D j = 2.5 and J = 13.3. A comparison of the predictions based on the corre- 
lations developed during this study with those of Reference 9 show somewhat 
closer agreement with the measured data using the techniques developed in 
this report. 

IV CONCLUSIONS 

A. Correlation Parameters 

The mixing efficiency and temperature distribution downstream 
from a row of multiple dilution orifices can be adequately predicted as a 
function of downstream distance over the range surveyed on this study provided 
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Conclusions (cont.) 


only that three independent variables are known: 

(1) The jet to primary stream momentum flux ratio, J 

(2) The nondimensional diluent orifice diameter, (H/Dj)"^'° 

(3) The nondimensional diluent orifice spacing, S/D^ 

This set of independent variables will allow predictions to be made for the 

following parameters: 

( 1 ) 

( 2 ) 

(3) 

(4) 

B. Model Precision 

The correlations developed during this study can be used over 
the ranges of variables given in Table I with reasonable confidence that the 
predictions will be within the one sigma standard error of prediction value 
given for each correlation in Table II. Extrapolation somewhat beyond the 
range of momentum flux ratios and downstream distances listed in Table I 
should yield reasonable predictions. However, extrapolations beyond the 
specified ranges of orifice size and spacing should be done with caution. 

That is, the correlations given will not reduce correctly to the limits 
of a slot jet or a single jet. Direct use of these correlations for combustor 
applications involves the implicit assumptions that the range of density ratios 
and turbulence levels surveyed during the test program of reference 1 were 
adequate to have allowed characterization if a significant influence existed. 


The mixing efficiency (energy exchange efficiency), Ey 
The jet velocity and temperature centerline penetration 

V°J> V D J 

The maximum nondimensional temperature values in the 
centerplane 

Shape factors which allow the entire temperature field 
to be predicted from the assumed Gaussian profile shape 
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TABLE I 

SUMMARY OF DATA RANGES 


PARAMETER 

NOMINAL 

RANGE 

Momentum Flux Ratio, J 

• 5.0 - 80.0 

Flow Rate Ratio, W,/W 

.04 - .60 

Density Ratio, p j/p to 

1.6 - 2.7 

Velocity Ratio, Vj/V w 

1.59 - 5.33 

Duct Height/Jet Diameter, H/Dj 

5-20 

Jet Spacing/Jet Diameter, S/Dj 

2.5 - 7.5 

Jet Spacing/Duct Height, S/H 

.125 - 1.0 

Downstream Distance/Duct Height, X/H 

.125 - 2.0 

Downstream Distance/Jet Diameter, X/Dj 

1.25 - 30. 

Primary Stream Reynolds Number 

.3 - .8 x 10 5 

Primary Stream Temperature 

450 - 750°K 

Primary Stream Velocity 

15 m/sec 

Jet Velocity 

25 - 121 m/sec 

Jet Temperature 

290° K 
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TABLE II 

SUMMARY OF CORRELATION EQUATIONS 

STANDARD 




PARAMETER 


Minus-Side Minimum Temperature 
Difference Ratio 


min, cent 

O' . 
c,cent 


Plus Side Half Width 


l/2,cent 

d j 


Minus Side Half Width 


Midplane to Centerplane 
Theta Ratio 


c.mid 


c,cent 


Off-Centerplane to Centerplane 
Theta Ratio 


c,cent 


Midplane to Centerplane 
Penetration Ratio 


c,mid 

^c,cent 


Off-Centerplane to Centerplane 
Penetration Ratio 


c,cent 


TABLE II (cont.) 


CORRELATION EQUATION 


[1.0 -e ] 

-.3 -1.4 .9 

c* = 1.57 (J) (S/Dj) (X/Dj) 


STANDARD 

ERROR OF PREDICTION 



.18 -.25 .5 .5 

0.162 (J) (S/Dj) ’ (H/Dj) (X/Dj) * 0.6598 

.15 .27 .5 .12 

•20 (J) (S/Dj) (H/Dj) (X/Dj) 0.5503 

-d 

[1.0 - e ] . 0.1120 




TABLE III 


MATRIX OF TEST CONFIGURATIONS AND MOMENTUM FLUX 
RATIOS USED TO DEVELOP CORRELATIONS 



2.5 

S/D 

3.75 — d- 

5.0 

7.5 

n/Dj 

5 

6-39^ 

N.T^ 

6-60 

N.T. 

7.5 

6.-60 

6-26^ 4 ^ 

N.T. 

N.T. 

10 

6-60 

6-30^ 

6-60 

6-60 

15 

1 4- 60 ^ 3 ^ 

N.T. 

6-60 

N.T. 

20 

6-60 

N.T. 

N.T. 

N.T. 


(1) No Tests Conducted with J Greater than 39 

(2) N.T. = Not Tested 

(3) J = 6 Test not Used - Invalid Thermocouple Data 

(4) J = 60 Test Not Used - Stored Test Data Could Not be Recovered 

Actual S/D. = 3.54; H/D. = 7.07 
J J 

(5) No Tests Conducted with J Greater than 30 
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TABLE IV 


MATRIX OF TEST CONFIGURATIONS AND MOMENTUM 
FLUX RATIOS USED TO ILLUSTRATE CENTERPLANE DIMENSIONLESS TEMPERATURE PROFILE! 


S/ Dj 


H/ Dj 

2.5. 

3.75 

5.0 

7.5 

5 

0=13 

N. t! 1 2 ^ 

0=27.2 

N.T. 

7.5 

Not ^ 
Used 

Not 

Used 

N.T. 

0=6 

N.T. 

10 

0=25 

0=25 

0=26 

0=60 

0=25 

15 

0=57 

N.T. 

0=60 

N.T. 

20 

0=25 

N.T. 

N.T. 

N.T. 


(1) Not Tested 

(2) Tested But Not Illustrated 
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TABLE V 


MATRIX OF TEST CONFIGURATION AND MOMENTUM FLUX RATIO 
USED TO ILLUSTRATE LATERAL PLANE DIMENSIONLESS TEMPERATURE PROFILES 


S/Dj 



2 , 5 

3.75 

5.0 

7.5 

H/ Dj 





5 

0=13.3 

0) 

Not 

N.T. 


X/H=l .0 

N.T. 

Used 


7.5 

Not ^ 

Not 

N.T. 

N.T. 


Used 

Used 



10 

0=25 

Not 

J=6 

J=25 


X/H=l .0 

Used 

J = 26 
0=60 
X/H=l .0 

X/H=l .0 

15 

0=57 

N.T. 

0=24.7 

N.T. 


X/H=l .0 


X/H=.25 


20 

J= 25 
X/H=l .0 

N.T. 

N.T. 

N.T. 


(1) Not Tested 

(2) Tested But Not Illustrated 
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FIGURE 1. MULTIPLE JET STUDY COORDINATE SYSTEM 













Velocity Penetration Depth, Y u /D 







Velocity Penetration Depth, Y /D 



FIGURE 5. COMPARISON OF EFFECT OF MOMENTUM FLUX RATIO ON 
JET VELOCITY PENETRATION AT X/D, = 10 




Y^Dj/C-539 J (S/Dj) (H/Dj) 




Heasured Centerline Temperature Difference Ratio, & 


FIGURE 7 „ JET CENTERPLANE CENTERLINE THETA CORRELATION 
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Downstream Distance, X/Dj 
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FIGURE 8. JET CEHTERPLANE PLUS SIDE 0 _ /9 CORRELATION 

1 /z ,cent 
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FIGURE 9. JET 'CENTER? LANE PLUS SIDE 0 ]/2 cent CORRELATION 
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FIGURE 11. 

CENTERPLANE TEMPERATURE PROFILE COMPARISONS 
J = 25.2, S/D. » 2.5, H/D. = 10.0 
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FIGURE 12. 


CENTERPLANE TEMPERATURE PROFILE COMPARISONS 
J = 26.8, S/D, = 5.T, H/D, = 10.2 
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CENTERPLANE TEMPERATURE PROFILE COMPARISONS 
J = 27.6, S/D, = 7.7, H/D, = 10.3 
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CENTERPLANE TEMPERATURE PROFILE COMPARISONS 
J = 61.9, S/D, = 5.1, H/D, = 10.3. 
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FIGURE 20. CENTERPLANE TEMPERATURE PROFIL 
J = 27.2, S/Dj -5.1, H/Dj « 5 
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FIGURE 23 LATERAL PLANE TEMPERATURE PROFILE COMPARISONS 

J ■ 61.9* S/D, = 5.1 » H/D, = 10.3, X/H = 1.0 
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FIGURE 26. 

LATERAL 
J = 25, 

PLANE TEMPERATURE PROFILE COMPARISONS 
S/D, - 2.5, H/Dj = 19.9, X/H » 1.0 
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LATERAL 
J * 13.3 


PLANE TEMPERATURE PROFILE COMPARISONS 
, S/Dj - 2.5, H/Dj • 4.9, X/H = 1.0 



APPENDIX 

TEMPERATURE FIELD PROGRAM 


l. 

2 * 

3 , 

9 . 

5 , 

P, 

7 , 

e. 

io, 

n, 

12 , 

* 13 , 
19 . 

15 , 

16 , 
17 , 

te. 

19 , 

20 , 
21 . 
22 . 
23 , 

. 29 , 

25 , 

26 , 

27 , 

28 , 

29 . 

30 . 

31 . 

32 . 
35 . 
39 , 

35 , 

36 , 

37 , 

38 , 

39 , 
90 , 
91 . 

92 , 

93 , 
« 9 . 
« 5 . 

96 . 

97 , 
<48 , 
99 , 

, 50 , 

51 , 

52 , 

53 , 
59 , 


C 

G THIS PROGRAM will USE THE EQUATIONS DEVELOPED DURING THE! MULTIPLE 1 
C JET STUDY, NAS 518 Q 26 , TO DEFINE THE: THERMAL) F I ELD DOWNSTREAM OF! 

C MULTIPLE JET INJECTION PORTS 

C 

C 

REAL NWIDTH 

COMMON / DIM / YCU 0 ),YV< 10 ),YH( 20 ),XH( 6 ),X 0 J( 10 )»ZS( 6 )i 
*nCAPUO>,TiMAX( 10 ),HAFPOSUO),HAFNEG(IO),PWlDTHC!Q) f ETUO)r 
*TWlDTH(lO),NWlOTH(IO)»TZ£( 6 , 2 n. YZZ( 6 , 2 l)iDEL( 10 )#Hr 

* THETA ( 10 . 21 , 20 ).XTHETA( 10 . 2 i. 20 ).TnL£U 3 ),YO( 20 )#YMID( 20 ),TMIDUO 
* ) , TQAR ( 1 0 ) 

COMMON / SINGLE / HD,RRHO,CD,RVEL,5DJ,RWD0T,HDJ,TmAL,XJ,S0, 
ftRTEMP 

COMMON / RPLOT / XPLOT( 22 ),YHPER( 22 ), I TAB,NDSTRM, NOPLOT ( 6 ) , YFIRST, 

aydel.xfirst.xdel 

C 

H=«,0 

c ***AXIAL DIST,, Y LOCATION# AND LATERPLt LOC AT ION ARRAYS *** 

DATA XH/,l£ 5 , , 250 , , 500 , 1 , 00 , 1 , 50 , 2 , 00 / 

DATA YH/, 039 # , 001 # , 127 , , 179 # . 221 #, 267 ,, 31 «#, 36 !#,« 08 # , 959 #, 501 # 
*, 598 # , 595 , , 691 , , 688 , , 735 , , 782 # , 828 , , 075 , , 922 / 

DATA ZS/.O, , 2 # , 9 , , 6 , , 8 , 1 , 0 / 

namelist/ input/hd, cd. xj»sd,sh,rrho', tideal’»rvel,rtemp,rwdot,IPRnt 
*, ITAH,NDSTRM,NQPL 0 T, YFIRST.YDEl, X F I RS T , X DEL , I PLOT 

1 FORMAT ( 1 3 Ab ) 

REA 0 ( 5 , 1 ) TITLE 
lF(CD.LT, f 0 I)COs ,62 
READ{ 5 , INPUT, END * 2000 ) 

WRITE( 6 , INPUT) 

c *** CONVERT FROM X/H AND Y/H TO X/OJ AND Y/DJ. 

2 FORMAT UHl , 15A6J 
DO 50 1=1,6 

50 XDJ(I) 5 XH(I)*H 0 /SQRT(CD) 

DO 60 1 * 1,20 

60 YDCI )sYHU)*HD/SORT(CO> 

IF (RWDOT ,EQ, 0 , ) WRITE ( 6 , 5 ) 

1 FCRWDUT ,EQ, 0 ,)RWDOT 9,20 
SDJ = St>/SQRT(CD) 

HDJ=HD/SURT (CO ) 

IFfRHHO.EQ.O, JRRH 0 * 2,2 
IF(SH t EQ, 0 .)SHsSO/HD 
IF (T IDEAL ,E 3 , 0 ,) TIDE AL»**DOT 
DO 100 1 = 1,6 
C *** ET EQUATION *** 

-N2 = “,681{5.(XDJ(I)**«99)*(SDJ***99)*(H0J**"1,0)*(XJ**,9J ) 

£T( I) = 100,0* (1 ,0-EXP (FN2) ) 

C *** CAP THETA AND MAX THETA EQUATIONS *** 

EXNsSQRT <SH/( 1 , ♦ SH ) } 

nCAP(l)s( 1 ,SJ6*UJ***,9 )*(XDJ(I)*a-1.15U**(£XN) 

T 1MAX(1)*1 1 C AP ( I ) M 1 ,*T IDEAL ) +T IDEAL! 

C *** PENETRATION EQUATIONS 

c *** THERMAL *** 

YC(I)=.559*(XJ** ,25) MSDJ**, 19)* ( HD J * * , 3 9 ) * ( XO J ( I ) * * • 1 7) 




FIELD 


55, 

56, 

57 , 

58, 

59, 

60, 
61, 
62, 

63, 

64, 

65, 
86 , 

67, 

68 , 

69, 

70 , 

71 , 

72 , 

73 , 

74 , 

75 , 

76 , 

77 , 

78 , 

79 , 
00. 
81. 
82 , 

83. 

84, 

05, 

86 , 

87 , 

88 , 

89, 

90, 

91, 

92 , 
95. 

94, 

95, 

96, 

97, 

98, 

99, 
100 , 
101 . 
102 , 

103, 

104, 

105, 

106, 

107, 

108, 
109, 

no, 

111 , 


EtX«EXP(*(XrtCl J **2 I 0)*(l,/SH»CCXJ«* t 5)/3,5l)/ll,0) 
YC(I)MC(I)*EEX 

c 

C *** VELOCITY *** 

VV(I)s,549ft(XJ**,12)*(SDJ**» 25) * (HOJ ** e 57 ) *(XDJ(!)**,IB) 

C 

C *** PLUS AND MINUS sioe half values 

c 

FX*,038*(XJ**1 ,62)* (X0J(I)**1 , l)*(H[>J**«3,67)*(8DJ*M ,5) 
HAFPOS(I)M ,«0,5*EXP(»FX) 

c 

FX»1 # 57*(xJ*a»,3J*(SDJ* < »»1(,4)*(XDJ(I)**,9) 

HAFNtGUJM ,-0,5*EXP<-FX) 

c 

c *** HALF WIDTHS 
C 

Tw|0THU)* t 3578*{XJ*A,17)*(HDJ*w,5)*(XDJU)*s3) 
PwIDTH(l) a ,l623*(XJ**,18)*(HDJ** t 5)*(XDJ(I)**,5)*(SDJA*»,25) f 
NWlDTH(n*TWlDTH(n-PrtIOTH(l) 

I F ( C YC ( I ) *P WIDTH { ! ) ) ,GT ,HOJ ) PWIOTH C I ) «HDJ*YC ( I ) 

IF((YC( D-NWIDTHCD) ,LT,0,) NWIDTH(1)bYC(I) 

100 CONTINUE 

c 

c *** END OF CENTER PLANE PARAMETER CALCULATIONS 
C 

C« *** BEGIN OFF CENTERPLANE EVALUATIONS • DATA INDICATES THE! RATIO' 
c= OF • ♦ AND • THETA HALF VALUES TO THETA MAX OFF CENTERPLANE IS. 

C ESENUAUY EUQAL TO THE CENTER PLANE RATIO, ALSO, THE OFF 

C CENTERPLANE HALFWIOTHS ARE EQUAL- TO THE: CENTERPLANE HALFHlDTHS 

C 
C 

c 

DO 150 1=1,6 

Fxa.4565 * (XJ*a ,528)* (SDJ»*- 1,529) MXDJU)**, 828) 

TM ID ( I ) M ,-EXPT-FX) 

FX3,227*CXJ**,67)*(SOJaa-1.0)*(XDJ(I)*«,S4) 

150 YHID(I)=t ,»EXP(-FX) 

C 

C. *** CALCULATE OFF CENTERPLANE! THETA ANO’ y MAX *** 

DO 200 1=1,6 
DO POO K=l, 6 

TZZ(I,K)M.M1*-TMID(I >)*Z$(K)**2 
20 0 YZZ(I,K)*1 # -U,-VMID(I))*2StK)A*2 
C 

C *** NOW wlU ASSUME GAUSSIAN DISTRIBUTION TO: GET FLOW FIELD **** 

C FIRST REPEAT THETA AND Y MAX VALUES FOR HALF SPAN FOUR TIMES 

PI) 250 1 = 1,6 
N = 0 

DO 210 *=7,1J 
Kl=5-N 

TZZC,K) = TZZU»Kt ) 

YZZII,K)=YZZ(I,KI) 

210 Ns N ♦ 1 
Nc2 

DD 220 K=12,16 
K 1 *N 

TZZU,K) = TZZU,K1) 




FIELD 


112 * * 
111 , 
116. 

115, 

116, 
117, 

lie, 

U9, 

120, 

121 , 

122, 

123, 

12a, 

125, 

1 2b , 

1 ?7 » 

1 2S , 

129, 

130, 

131, 

132, 
113, 
13a, 

115, 

116, 
117, 
136, 
139, 

iao, 
ui , 
162 . 

165. 

iaa, 

H5, 

tab, 

167, 

lae, 

169, 

150, 

151, 

152, 

155, 

156 , 

155, 

156, 

157, 

158, 

159, 

160, 
161 , 
162, 
163, 

166 , 

165, 

166, 
167, 
166 , 


Y2ZUpKmZZ(I»K|) 

220 N*NM 

N*0 t 

DU 230 K*J7,2l 
K 1 sSc'N 

T22(l»K)sTZZ(IfKl) 

YZZ<I,KJ*rZZ<l»K|) 

230 N*NM 
250 CDNTINUEi 

c 

c 

c 

c 

e 

DO 500 M*t,6 
DO 500 Ksi»2i 
DO 500 1*1,20 
YMAXbYZZ(M,K}*YC(M) 

IF ( YD ( I ) ,GE , YMAX ) GO TO 600 
Y I 8YMAK-YD C I ) 

TMIN=2,*HAFNEG(M)*T1MAX(M)*T1MAX(H) 

XEXPbEXP(.A\.ODC2,)*(YI/NW1DTH(M>)**2) 

THETA(H,K,l)PtZZ(HrK)*((TlMAX(H)-THlN)*XEXP+TMlN) 

GO TO 500 

aoo yi sydu j-vmax 

TMINs2,*HAFP0S<M)*T1MAX(M)»T1MAX(M) 

XEXPsEXPt-ALaGia.^m/PWlDTHCM))**?) 

THETA(M,K f I )STZZ(M,K)*{(T1MAX{M)-TMIN) *XEXP4TMIN) 

• 5Q0 CONTINUE 

c 

c THE FLOW FIELD HAS BEEN DEVELOPED OYER A TWO 8 SPAN FROM CENTER: 

C PLANE OF ORIFICE TO CENTER PLANE OF ORIFICE, NOW TRANSPOSE' 

C TO A FLOW FIELD THAT GOES PROS MIDPLANE: TO MIDPLANE OVER 25- SPAN 

C " 

C 

C : 

DO 700 Ms 1 , 6 
DO 700 Ksl, 16 
DO 700 IM,20 

700 XmTA<M,K,I}sTHETA(M,K*5,I) 

DO 720 *M,6 
DO 720 K=l7,2l 
DO 720 1*1,20 

720 XTHETA(M,K f I )tTHETA(M,K«l5, 1 ) 

C *** GET average THETA 
DO 750 Msl, 6 
T 6 AR ( M ) *0 , 

DO 750 KM, 21 
DO 750 1*1,20 

750 TBAR(H)3T6AR(M)+XTHETA(M,K,I)/20,/2l # 

C 

C PATTERN FACTOR 

C 

DO 760 MM , 6 

760 DEL(H)®T(iAR(M)/U,-TBAR(H)) 

C 

fFUPRNT.EO.OjCALL PRINT1 




FIELD 


169, 2000 CONTINUE 

1 70 s 3 F0RMAK1OX, '**« FLOW RATIO' INPUT AS ZERO# SET EQUAL TO 0,20 

t. 7 I , C 

172, C 

17 3* C 

179, C‘ 

175, IFUPL0T,EQ,0)CALL PLOT! 

176, END 

END ELT, TIME! O,«160 SECONDS, 
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PLOT i 


1. 

2. 

$ , 
«■ 

5 , 

6 . 
7 < 
6, 
9 0 

10 . 

11. 

12, 

13, 

14, 

15 , 

16 , 
17, 
19. 

19. 

20, 
21. 
22. 
23. 
2“. 

25, 

26, 
27, 
26. 

29, 

30, 

31 , 

32, 

33, 
3«. 

35, 

36, 

37, 
36. 
39, 
«Q. 

91, 

92, 
43, 
9 4, 

45. 

46. 

47. 
46, 
«9, 
50, 
51 , 

52. 

53, 
5«, 


SU9RDUUNEi PLOtl 
REAL NfUDTM 

COMMON / DIM / YCUO),YniO),YH(20>,XHC6),XDJUO),Z$C6),71CAP(lO>, 
*TlMAXnO),HAFPOS(lO),HATNEC(lO),PWIOTH(lO),eTUO)#Twn>TH(10). 

*N« IDTH< 10 ), TZZ (6, 2 1 ), YZZ (6, 2 1), DEL C 10 ),H, THETA (10,21,20)* 

• X7HETAU0, 21,20 ), TITLE (13), YO( 20), YMID(20),TMlO(2O),T8ARUO) 
common / single / hd,rrho,cd,rvel,$dJ/Rwdqt,hdj,tideal,x«J,sd,rtemp 
COMMON / RPlOT ( XPL0T(22),YHPER(22),ir A9, ND5TRH, NOPLOT (6),YF|RST, 

*YDEL,XFlRST,XOEL 

1F(NOSTRH,Eq,0)NDSTRHs« 

CALL PLOT S ( 0 , 0, 7 ) 

CALL PLOT ( 7 , # 1 . 0 * ®3 ) 

e *** (nm)c number of point locations in duct height *** 

C' *** (NDS TftM) = DOWNSTREAM LOCATION OP LATERAL PLOTS *** 

C »** (NOPLOT)s DOWNSTREAM OR LATERAL- LOCATIONS TO BE DELETED *** 

NM 3 2 0 

IF(ITA3,EQ.1)CALL SYMBOL, 10, 9, 7,0, ,20# ' CENT ERP LANE TEMPERATURE PROF 

* I LE COMPARISONS* ,0,, 45) 

IF(I7AB.EQ,2)CALL SYMBOL (0,1# 7,0, ,20# 'LATERAL PLANE. TEMPERATURE! PR 
•QFILE COMPARISONS' ,0,, PS) 

CALL STweOLl 3,8#6,5. .10, * J*' »0,,2) 

CALL SYMBOL(5 t 7#6,5,,10#'52DJ»',O,#5) 

CALL SYMBOL (8, 1,6 ,5# , 10, ' H/D J * ' , 0 , , 5) 

CALL NUMBER(U,4,6.5,,10,XJ,0,,1) 

CALL NUMBER (6, 7, 6,5, , 1 0 , SO J , 0 . , 2 ) 

CALL NUMBERS, 1,6, 5,, 10, HDJ,0,,i) 

IF(ITA8,EQ,2)CALL 5YMBOL(9,8,6,0, , 10, ' X/H* ♦ , 0 , , 4 ) 

IF(ITAB,EQ,2)CALL SYMBOL ( 7, 2, 6.0, , 10, ' X/O J» » , 0 , , 5 ) 
IF(iTA3,fcQ„?)CALL NUMBER (S, 7,6,0, , 1 0 , XK ( NOS T RM ) , 0 , , 3 ) 
IF(J7AB,EQ„2)CALl NUMBER (8 , 2, 6, 0, , 1 0, XDJ (NOSTRM) , 0, , I) 

CALL AXIS 10. ,0.# ’PERCENT OF DUCT HE I GHT » , 22, 5 , , 90 , , YF IR ST , YOEli) 

YHPER(NMM JsYFJRST 

YHPER(NM+2)sYDEL 

XPLOT (NH+l )*XF 1RST 

XPLOT (NM*2)sXDEL 

MMsO 

NMQs i 

10 DO 15 Mb 1 f 6 

IFtNOPLOT (NMD),tQ,M)GO TO 15 
DO 13 1*1 ,20 
YHPER(I)5100,*YH(I) 

IF(ITAB,EQ, 1) XPLOT (I)«X THETA (M, 6,1) 

IF C|TA0,EO,2)XPLOT (I) = XfHETA(VI0STRM,M*5>I) 

13 CON ? I nUE 

IF (MM.EU, 0)GO TO 14 

CALL AXIS(0,,0.#» »,5,S,,90.,YFIRST,YDEt) 

14 CALL AXIS(0,,0., »THETA»,'.5,3,,0,,XF1RST,XDEL) 

MM=MMfl 

CALL LINE (XPLOT , YHPER,NM, 1,0,1) 

IF UTA3.ED, HCALL SYMBOL < 0 , 7 , 5, 5 , , l 0 , ' X/DJ» ' , 0 . , 5) 

IF( I TA3.EQ, I) CALL NUMBER (1.7 ,5,5, ,10,XDJ(M),0.,1) 
lF(ITAa,ea,J)CALL SYM8QL(0.7,5 f 2,.iO,»K/H=*,O,,4) 

IF (I TAB, £ Q t DC ALL NUMBER (1 ,6,5.2, , 1 0, XH< M) , 0, , 3) 

IF ( ) T A B . E 1 , ?) r A| L SVM»0U to. 0,5,5, . 10, ’?;/$=' ,0,,5) 
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plot* 


END EL T , 


5 5, 

56, 

57, 
56 , 
5 P„ 
60, 
M, 
62, 
65 , 
6 R, 


IFUTAa t eo 0 2)C*LL NUMBER (1 .9,5,5, ,10, 25(H) ,0, ff 1) 

IF ( ITAB.ep.a, ANO,M,EQ, I ) CALL S YMBOU ( 0 , 6 , 5 , 3# , 10, * (CENTERPLANE) ! #0, 
*t 13) 

IFUTA8,EQ,2,ANO,M t E0,6KAU 3YM80U ( 0 , 0 , 5, 3; # 1 0 , » (MIDPt,ANE) 1 > 0 , « l Q 
* ) 

CALL PLOT (3,5, 0,0, «3) 

15 If (NOPLOT < NMD ),tQ,M>NMO»NMO+ I 
CALL PLOT (l<l,»*i,»999) 

RETURN 

END 


TIME* 0,UM SECONDS, 


- 65 - 



PRINT I 


1 » 

2. 

3, 

6, 

7, 

&, 

9, 

10 , 

n, 

12 , 

13, 
I®. 

15 , 

15 . 
17« 

16. 
1 ’. 
20 , 
21, 
22. 
23. 
20. 

25, 

26, 
27, 
26 , 

29. 

30. 
3 1 o 

32. 

33, 
30, 

35, 

36, 

37, 

36, 

39, 

00. 

01 . 

02. 

03, 

00 , 

05 , 
0 6 . 
07, 

06. 
09. 

50. 

51. 

52 i 

53, 

59 . 


SUBROUTINE PRINTS 
REAS, NdlDTM 

COMMON / DIM / YC(10),YVUO),YH(20),XH(6),XOJ<tO),Z$(b), 
*UCAPUO), TlMAX(lO), HAFPOSCIO), HAFNEGUO), PWlDTH(lO), ET(IO), 
*TWJD7H( 10),N«IDTH(10), 7ZZC6.21 ) , YZ2 (b, 2! ) , DEL ( 1 0 ) , H, 

*THET All 0,21, 20), XTHETA (10, 21, 20), TITLE (13), YD(2G),YMID(20),TMJD(20 
*),T9ARU0) 

COMMON / SINGLE / HO . RR HO , CO , R VEL , SDJ , RwDQT , HDU , ? 1 Q£ AL , X J , SO , 
APUMP 

DIMENSION ZZS(2i ), YS(20) 

HSlaH^a.Sa 

DATA 7^5/ Op,, 1,, 2,, 3, ,4,, 5#, bp g7,pQf, 9^1,0, 1,1,1, 2, 1,3, 1,4, 1.5, 

M ,6, 1 ,7, 1 „0, 1 ,9,2,0/ 

WRITE <6, 17) 

w*ms:(6, unme 

WRI TE(b, 1 ) 

WRITE(6,2)XJ,$D 
WR1 TE{6,3)RTEMP,H0 
WRITE(6,9)RRHU,CD 
WWI TE ( b,5)RVEL»S0J 
WRIT£(6|,6)RwD0T, HDU 
WRITE(b f 7)TI0EAL#H3I.H 
WRI T£ (6,0) 

WRIT£( 6 , 9 ) 

WRITE(6 # 20) 

DO 100 M=t,6 

100 WRITE (6, lO)xH<M},XDJ{M),ET<M),TiCAP{H),71MAX(M),YC(M),YV(M),T&Aft(M 
*),OEL<M) 

DO UO 1*1,20 
110 YS(I)=YH(I)*H0 
00 500 
WRITE(b,m 

WRIT£(6, 12>XH(M) ,XDJ(M) 

WRm(6,13)<ZZ$(K),Ksl,lU 
WRITEC6, SU) 

00 150 1*1,20 

ISO WRITE ( 6 , 15)YH(n,Y0(I),YS(I),(XTHETA(M,K#I),K*l, U) 

WRI T£(b, 12)XH(M),XDJ(M) 

WRITE (b, U) (ZZS(K) ,K*U.2l) 

WRIU(b, 14) 

00 160 1*1,20 

160 WRI IE (6, 15)YH(nrYDU),YSCn,(XTHETA(H,K.n,KBU,2i| 

500 CONTINUE 

1 FORMAT (///20X, ****** OPERATING CONDITIONS *aa*a», T bQ , * **** DESIGN: 
*C0NDITI0NS »*»**'//) 

2 FORMAT (20x p »MOMENTUM FLUX RATIO* • , Fb, 2, T60, * OR IF ICE; SPACING, $/D« 
a »,T90,F6.3) 

3 FORMAT (20X, * TEMPERATURE RATIO * » , F6,2, 760, » OR ZMCE- SIZE # H/D. * 
a 1 #T90#F6,3) 

4 FORMAT (20x, 'OENSITY RATIO * *,£6.2,760, 1 ORIFICE- DISCHARGE! CO£i 

*F8»T90fFb t 3) 

5 FO«MAT(20X, ^VELOCITY RATIO * ‘ , F6.2, T60, * EFF£CU V£ SPACING, S/D 

Aj* , T90,F6,3) 

6 FURMAH 20 Y,'FLOW RATE RATIO* • , F b, 2 , T 60 , ' FFF FCT I Vf ORIFICE 312 
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PR l NT 1 


55, 

56, 

57, 
56, 

59, 

60 , 
61, 
62, 
63, 
b«, 

65, 

66 , 
67. 
66 , 
69, 
7.0, 
7 t « 
72, 
7i t 
7«, 

75, 

76, 


*E, H/OJs *,T 90 ,F 6 ,J) 

7 FQRMAH 20 X, 'IDEAL THETA a * ,F 6 , 2 ,T 60 , »OUCT HEIGHT® *,T 90 , 

*F fa , 3 # 1 CM 'Fb. 3 ,‘ IN M 

e FORMAT W// 40 X ,'*****» MIXING AND CEMTERPLANE OAT A ****»**///) 

9 FORMA T ( 1 3 X , ' DISTANCE', T 27 ,»OISfANC£iST« 0 , ‘MIXING EFF « , T 52 , ' CAP. TH 
*ETA' ,T 64 , ’MAX THETA' ,T 76 , 'PENETRATION' ,TS 8 , 'PENETRATION ’, 1100 # ‘ AYE 

* THETA*, T 112 , ‘PATTERN* ) 

20 FORMATU 7 X,'X/H*,T 30 ,‘X/DJ*,T 45 ,*ET*,T 75 ,* X/DJ ( TEMP) «, TBS » ’ X/D 
*J (VEL) ' , T 115 , 'FACTOR ' ) 

10 Fl)RMATU 0 X, 9 t 2 X,F 10 .R)) 

U FORMAT ( JHI ,-T 50 r ' *** table OF THETA VALUES A**'//) 

12 FORMAT (IH 0 »TU, f X/H' ,F 6 ,«,T 30 , ‘X/ 0 J*#PB, 4 /) 

13 PORMATtTR, *Z/$s*,T 20 , 11 F 10 ,«) 

14 FORMAT! T 5, ' V/H»,Tll, * Y/O J * , T 1 7 , ' Y/S » ) 

15 FORHAT(F 5 f 3 ,F 0 t 3 ,F 6 , 3 ,T 2 O, UFiO,«) 

16 FORMAT 120 X, ' ****** ‘ , 1 3 A 6 , T 11 0 , ‘ ******* ) 

17 F OR MAY{lHi,////// 2 QX» *********** AEROJET LIQUID ROCKET COMPANY HUt 

* T l PLE SET INJECTION FLO* FIELD PROGRAM *********** /// 45 X, * ( 5 EVEL 
*OPEO ON NASA LEWIS CONTRACT NAS 3 - 19026 )'///) 

RETURN 

END 


3X0T FIELD 


SInPUT 

HD 

8 

,4000000 OE +01 

CD 

8 

,661 30000E+00 

XJ 

S 

, 1 5320000E+02 

SD 

8 

,20000000E+01 

SH 

S 

, 50000 0OOE + OO 

RRK3 

8 

, 22 1 00000F + 0 1 

T I DEAL 

■ 

,27590000 E +00 

«V£L ! 

8 

, 20500000 E + 0 1 

RTEMP 

8 

.22100000E+01 

RWOOT 

* 

, J500000QE+00 

IPRMT 

3 

+ 0 

I T A 9 

8 

+ 1 

NDSTRM 

8 

+ a 

noplot 

S 

♦ 1 , 

Y FIRST 

- 

*C, 

,ooooooooe+oo 

Ydel 1 

8 

.20000000E+02 

XFIRST 

8 

,OOOOOOOQE+O0 

XDEL: 

8 

.2S00C000E+00 

i plot 

= 

+ 0 

$£ NO 




SAMPLE .INPUT 


♦ 5, 

*0 


♦ 0 , 


*Qf 
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